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ELASTIC CONSTANTS FOR BENDING AND TWISTING OF 


CORRUGATION- STIFFENED PANELS* 
By W. Jefferson Stroud 


SUMMARY 


Theoretical formulas for five orthotropic elastic constants associated with 
pure Lending and pure twist are presented for a corrugation- stiffened panel. The 
constants are the bending stiffnesses parallel to and normal to the corrugation 
direction, the twisting stiffness, and the Poisson's ratios associated with 
bending. The formulas are given in terms of the dimensions of the corrugation 
pattern for the general case of different corrugation and cover- plate materials. 
Bending stiffnesses measured on panels with three different corrugation patterns, 
and twisting stiffnesses measured on panels with two different patterns, compare 
favorably with the theoretical results. 


INTRODUCTION 


A type of structural panel used in the skin of the X -15 and having potential 
use for the skin of hypersonic aircraft and reentry vehicles consists of a corru- 
gated metal plate attached to one side of a flat metal plate. Such a panel may 
be thought of as a plate stiffened by corrugations and is referred to in this 
report as a corrugation- stiffened panel. This type of stiffened panel is easy to 
fabricate and for those applications where the primary load- carrying function 
requires a high bending stiffness in one direction, it provides an efficient use 
of the material. Its main disadvantage is that the bending stiffness in the 
other direction is extremely small. Thus it may be susceptible, in some applica- 
tions, to undesirable effects such as panel flutter. 

In order to predict the structural characteristics of a corrugation- 
stiffened panel (its panel- flutter characteristics as well as its load- carrying 
ability), information is needed about the elastic constants of the panel. These 
constants define the resistance of a loaded panel to deformation. As discussed 


*Part of the information presented herein was included in a thesis entitled 
"The Bending and Twisting Stiffnesses of Corrugation Stiffened Panels" submitted 
in partial fulfillment of the requirements for the degree of Master of Science 
in Engineering Mechanics, Virginia Polytechnic Institute, Blacksburg, Virginia, 
June 1962. 



in reference 1 , in which corrugated- core sandwich panels are considered, the con- 
stants include two transverse shear stiffnesses Dq^ and Dq_^, two tending stiff- 
nesses D x and Dy, a twisting stiffness Dxy, two stretching moduli E x and 
Ey, a shearing modulus G^, two Poisson's ratios q x and |iy associated with 
bending, and two Poisson's ratios n' x and n' v associated with stretching. 


In the present paper the five elastic constants associated with pure bending 
and pure twist, D x , Dy, D X y, (j. x , and p.y, are considered. Formulas are given 

for these five constants in terms of the dimensions of the corrugation pattern 
for the general case of different corrugation and cover- plate materials. Theo- 
retical values of D x , Dy, and D-^y are compared with experimental values for 

three corrugation patterns. For each specimen tested, the cover- plate material 
was the same as the corrugation material. 


SYMBOLS 


area enclosed by corrugation cell 


a ,b,d,e-j_, 

n.Pjr-L,^, 

y,e 

c l> c 2 } c 3 
D 

D X ,Dy 


D 


’xy 


D* 


V 


E 

G 

H 

I 


dimensions of corrugation pattern (see fig. 2) 

constants in solution of differential equation in appendix D 

Et 5 


plate stiffness per unit width, 


12 


i 1 - » 2 ) 


, in- lb 


bending stiffnesses per unit width of a beam cut from orthotropic 
panel in x- and y-directions, respectively, in- lb 

twisting stiffness of unit-width and unit-length element cut from 
orthotropic panel, with edges parallel to x- and y-axes, in- lb 

beam bending stiffness in y-direction per unit width for segment 
of panel of length 2b (defined in eq. (B19)) 

beam bending stiffness in y-direction per unit width for segment 
of panel of length 2a (defined in eq. (B27)) 

Young ' s modulus, lb/sq in. 

shear modulus, lb/sq in. 

quantity associated with shape of corrugation (defined by eq. (B12)) 

moment of inertia per unit width, in. 5 


2 



I 

J 

K 

M 

Mx> M y 


M- 


L xy 


n 

P 


area moment of inertia per unit width of cross section, in. 5 
torsion constant per unit width of cross section, m.-' 
curvature in x-direction 
moment per unit width, in-lh/in. 

resultant bending- moment intensities in x- and y-directions, 
respectively, in- lb/in. 

resultant twisting-moment intensity with regard to x- and 
y-directions, in-lb/in. 

integer 

in-plane load per unit width acting on corrugation and cover plate 
at weld line, lb/in. (see figs. 13 and 15 ) 

radius of curvature in y- direction, in. 


s coordinate measured along center line of sheet making up corrugation 

cell, in. 

t thickness, in. 

U strain energy due to bending, in- lb 

u,v,w r displacements in x-, and radial directions of middle surface of 

a cylinder, respectively 

functions arising from integrations (see eqs. (D9) and (DIO), 
respectively) 


w displacement in z-direction, in. 

X quantity associated with shape of corrugation (defined by eq.. (Bl4)) 

x coordinate, measured parallel to corrugation direction, in. 

y coordinate, measured parallel to cover plate and perpendicular to 

corrugation direction, in. 

z coordinate, measured perpendicular to cover plate, in. 

z distance from outer surface of cover plate to neutral axis of 

segment of length 2 a, in. 

7 X 0 shear strain, in. /in. 

5 extension of elements 


5 



e strain, in. /in. 

0 angular coordinate for elements BC and DEj from B to C for ele- 

ment BC, from D to E for element DE, radians (see fig. 13) 

H Poisson's ratio of material 




^xl 


dx2 


1 

2 


0 


Poisson's ratios associated with curvature in y-direction caused 
by M x loading and with curvature in x-direction caused by My 

loading, respectively 

Poisson's ratio in x-direction for segment of panel of length 2b 
(defined by eq. (D52)) 

Poisson's ratio in x-direction for segment of panel of length 2a 
(defined by eq. (D55)) 

distance measured along element CD, in. (see figs. 13 and 15) 

quantity associated with material and shape of corrugation 
(defined in eq. (D^3)) 

stress, lb/sq in. 

angular coordinate measured clockwise from the vertical, radians 
(see fig. 15 ) 

quantities associated with shape of corrugation (defined in 
eqs. (D4l) and (D^9)) 


Subscripts : 

c refers to corrugation 

exp experimental 

p refers to plate to which corrugation is welded 

th theoretical 

x refers to x-direction 

y refers to y-direction 

£ refers to 1-direction 

0 refers to circumferential direction 

1,2,. . . 5 refers to extensions in y-direction of elements AB,BC,. . . EF, 

respectively 



THEORETICAL RESULTS 


The bending and twisting elastic constants used herein are defined in the 
following pure-bending and pure-twisting relationships (see, for example, 
ref. 2): 


D. 


■V~r 


'y /d^w , d^w 

* — + Px — 


PxM-yydy^ dx? 


M- 


xy 


= D 


xy 


d^w 
dx dy 


The directions for M x , My, and Mxy - 
and, consequently, the directions for 
D x , Dy, and D X y - are indicated in 

figure 1. The stiffnesses Dx, Dy, 
and D X y are associated with the 
resultant moments M x , My, and M X y, 
respectively. When the transverse shear 
stiffness is assumed infinite, the dif- 
ferential equation for lateral deflec- 
tions of an orthotropic panel consistent 
with these definitions of the elastic 
constants is 



Figure 1.- Typical corrugation-stiffened 
panel. 


D x d\r + o/p t ^y°x \ d^w + D y d^w 
" ^xHy dx^ - \ ^ “ ^xH-jy dx^dy^ 1 “ H-xPy dy^ 


. . „ o'Ht , tit d^ 

^ Nx _ +N y_^ X y_ 


+ 2N-. 


where 


q lateral load, lb/sq in. 

N x ,Ny intensities of resultant normal force acting in x- and y-directions, 
respectively, lb/in. 

Nxy intensity of resultant shearing force in xy- plane, lb/in. 


Derivations of theoretical expressions for the three stiffnesses D x , Dy, 
and D X y and the Poisson's ratios p x and |iy for cross-section shapes of the 

type shown in figure 2 are outlined in appendixes A, B, C, and D for the general 
case in which the cover- plate material and corrugation material are different. 

For easy reference, results for the case in which the cover- sheet and corrugation 
material are the same are given by the following expressions: 
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For the bending stiffness in the x-direction: 

D x = El 

where E is the Young's modulus of the panel material, 

2 


( 1 ) 


T-§ 


+ (a+bHp.^i^ + (e 1+ r^ ^(a+i)t c + t c r^^) 

+ 2t c eiri 2 sin 0 + ei 2 rp0t c + ^ k3sin 2 0 + cos 2 0 + ^tp + ^- + j + ^ sin 0 - kt c 

+ t c r 2^(^| + ? " 1 ' \ 29 ) + 2tc e 2 r 2 2sin 9 + e 2 2r 2 9 t c + + (e 2 + r 2 ) 2 t c g 


( 2 ) 


and 


y = 


^-(b+a) + (t p + ( a+i ) t c + (t p + ^ + ri - - riS ^ n - 9 ) ri 0t c + (t p + ^ + j + kt c 


+ (tp + ^ + d-r 2 + r ^ n e j r 2 0t c + (t p t^-+dj 


gtc 


t p (b+a) + (a+i)t c + r]_0t c +kt c + r 2 0t c + gt c 


(3) 


All the symbols on the right-hand sides of equations ( 2 ) and ( 3 ) are dimensions 


shown in figure 2. 

In particular, note from figure 2 

that 


II 

1 

''V 

H 

+ 

c+ 

+ 

jo jo* 

(M 

, = d - ( r i + r 2 + e x ) 

(5) 

j = r-^1 - cos 0) 

(6) 

= d - ^(l - cos 0) 

(7) 



Figure 2.- Dimensions of corrugation 
cell. 


Throughout this report, the subscript c refers to the corrugation and the 
subscript p refers to the .plate to which the corrugation is attached. 
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For the bending stiffness in the y-direction: 


Dy = a + b 

V V 


( 8 ) 


where 


D y' - ^ x + Dp )( 1 ” 


„ _ E fe + M 3 , 


12 


CR— —Si—fl Li 4. 

>(l - ^2) l 1 “ '**) 

In equations (9) and (10), m- is "the Poisson's ratio of the panel material, 

Et c 5 

i 1 - 


(9) 

(10) 


B c = 


12 


Bp = 


Et p 5 


12 


jrr^j 


(id 

( 12 ) 


and 


where 


H = = 


X = i + ^r^ + rgj 9 + k + g - Hjr 1 2 (e - sin 0) + jk + — sin 0 
+ rg^nO - rg0 cos 0 + r2 sin 9^ + gdj 

2 

r x 2 (e - sin 0) + jk + sin 9 + rg^nO - rg0 cos 0 + rg sin 0^ + gd 


(13) 


r l 


3^30. _ 2 sin 0 + sln 29 ^ | + j 2 k + j^ 2 s i n 9 + ~ sin 2 0 + r 2 n 2 0 - 2nr 2 2 0 cos 0] 


+ 2m*2 2 sin 9 + (30 - 3 sin 9 cos 0-20 sin 2 0) + gd 2 


(14) 


The quantities p. x and p.y are defined subsequently but the product n x My I s 
usually small compared with unity and can be neglected in the calculation of Dy' 


and Dy" . 


L 



For the twisting stiffness: 



(15) 


where G is the shear modulus of the panel material and 


| |'b^t c + tpj + d|b + g - i + (r 2 - r-^ tan | +2^2 - r 2 2 jtan | +-(r 2 2 - r-,2 

j = ^ £ 

t5 + ij[ 1+6 ( r l +r 2) +IS+g ] 

For Poisson’s ratio associated with curvature in the y-direction caused by 
M x loading : 



p x = 


bp x i + ap 
b + a 


where 


and 


^xl - 


-HE 


D p 

D? X + " 


(17) 


( 18 ) 


Z = p ^i ^ + r^j + r^e^. sin 9 + r^sin q cos 9 + co 7 . P ^ 


e^ + r-j_ cos 9)' 


- ^e 2 + r 2 cos 9^ 


^ ^ - r 2 e 2 sin 9 - r 2 2 sin 9 cos 9 - g^e 2 + r 2 j - b^y - 


(19) 


For Poisson's ratio associated with curvature in the x-direction caused by 
My loading, the following relationship, which can be derived from the theorem of 
reciprocity (ref. 3 )* I s used: 

Py = Hx jr (20) 

u yL 

In the derivation of equations (l) and (l5)> the panel was considered to 
be made up of a series of cells or thin- walled tubes. In addition, for equa- 
tion ( 15 ), the angle of twist of each corrugation element was considered to be 
equal to the angle of twist of the panel as a whole. 
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EXPERIMENTAL VERIFICATION OF THEORY 


General 


The values of D 


and D were determined experimentally for three 


x' y’ ' xy 

corrugation patterns — the small vee ; the large vee, and the square — shown in 
figure 5. A total of nine panels were tested. Each panel had a nominal depth 
of l/2 inch. The small -vee corrugation had a nominal pitch of l/2 inch; the 

large -vee and the square corrugations 



(a) 


1.0 



(b) 


1.0 


.14 


.22 


.14" 


.5" 


Weld 

lines 


6 = 90" 


(C) 


(a) Small vee; = r 2 = 0.05". 

(b) Large vee; r i = rg = O.lV. 
(c) Square; r^ = r 2 = 0". 

Figure 3.- Corrugation shapes. 


each had a nominal pitch of 1 inch. 

The variation in depth and pitch for a 
specimen did not exceed ±4 percent 
except for the edge corrugations. 

Some of the panels were fabricated of 
AISI type 347 stainless steel; others 
were fabricated of Inconel X. Spot 
welds of approximately 0.2-inch spacing 
were used to join the corrugated plate 
to the cover plate. Both the small-vee 
and the large-vee corrugations had sin- 
gle weld lines; the square corrugation 
had a double weld line as indicated in 
figure 3* The materials and average 
dimensions of all the test panels are 
given in table I. (The values of E 
and |i for type 347 stainless steel 
were obtained from ref. 4. The values 


of, E, G, and p for 
obtained from ref . 5 . ) 


Inconel X were 


In order to assess the importance 
of Poisson's ratio in the determina- 
tion of the measured bending stiff- 
ness, Poisson's ratios li and u 

x y 

were calculated for the test panels by 
using equations ( 17 ) and (20) . These 
values are tabulated in table II. 

is less 


^y 


Note that the product 

than 3 X 10 “5 for all the test panels. 
Since small values of this product are 
typical of corrugation-stiffened 
panels, elementary beam theory can be 
used in the calculation of stiffnesses 
from measured deformations of the panel. 
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TABLE I.- AVERAGE DIMENSIONS OF EACH TEST PANEL 


Test 

panel 

Corrugation 

pattern 

P, 

in. 

a, 

in. 

in. 

gj 

in. 

in. 

a, 

in. 

r l> 

in. 

r 2 > 

in. 

0, 

deg 

"kpj 

in. 

to 

in. 

y, 

in. 



AISI type 

3^+7 stainless 

steel: 

•E = 28.6 X 10 6 ; p = 

0.30 





Small vee 

0.5290 

0 

0 

0 

0.4512 

0.4900 

0.0543 

0.0543 

68.79 

0.0150 

0.0150 

0.1862 

Dx 

Large vee 

1.000 

0 

0 

0 

.4598 

.4700 

.1380 

.1380 

52.18 

.0150 

.0150 

.1543 


Square 

1.000 

.1120 

.1390 

.2^90 

.4830 

.4830 

0 

0 

90.00 

.0150 

.0150 

.1772 



Inconel 

X: E = 

31.0 x 10 6 ; 

= 11.0 x 10 6 

p = 0 

29 





Small vee 

0.5000 

0 

0 

0 

0 . 4601 

0.5040 

0.0520 

0.0520 

70.73 

0.01064 

0.01064 

0.1900 

D y 

Large vee 

1.000 

0 

0 

0 

.4474 

.4750 

.1450 

.1450 

53.26 

. 01068 

. 01065 

.1514 


Square 

• 9950 

.1150 

• 1330 

.2488 

.5063 

.5063 

0 

0 

90.0 

. 01067 

. 01054 

.1812 


Small vee 

0.5020 

0 

0 

0 

0.4555 

0.4981 

0.0520 

0.0520 

70.37 

0.0106 

0.0106 

0.1869 

D xy 

Large vee 

1.006 

0 

0 

0 

.4507 

.4750 

.1450 

.1450 

52-95 

.0106 

.0106 

.1513 


Square 

.9960 

.1060 

• 1375 

.2545 

.4960 

.4960 

0 

0 

90.00 

.0106 

.0106 

.1795 


TABLE II.- CALCULATED VALUES OF POISSON'S RATIOS FOR EACH TEST PANEL 


Test panel 

Corrugation pattern 

^x 

^y 

q x ^y 


Small vee 

0:141 

7.70 X 10" 5 

1.08 x 10" 5 

D x 

Large vee 

.156 

15.2 

2.38 


Square 

• 175 

8 . 4 i 

1.47 


Small vee 

0.130 

3.16 x io "5 

0.413 x 10 " 5 

D y 

Large vee 

.147 

7.00 

1.02 


Square 

.169 

3.65 

.616 


Small vee 

0.131 

3.30 x 10‘5 

0.434 x io -5 

D xy 

Large vee 

.147 

7-04 

1.04 


Square 

.168 

3-77 

.632 


Experimental Verification of Formula for D x 

Specimens . - The three specimens — one with each of three corrugation 
patterns — used for the evaluation of D x were beams approximately 44 inches 
long and approximately 6 inches wide with corrugations parallel to the long 
dimension. Each specimen was fabricated of type Jk 7 stainless-steel sheet 
which had an average thickness of 0.0150 inch with a variation in thickness of 
±2 percent , 

Test apparatus and procedure .- The test setup with pertinent dimensions is 
shown schematically in figure 4. A photograph of the apparatus with specimen 
is shown in figure 5* The beam was supported on two knife edges and loaded at 
the ends so as to obtain a region of pure bending moment between supports. The 
supports were 28.47 inches apart, and loads were applied with lead shot in 
increments of. about 5 pounds to a maximum of about 60 pounds. Deflections of 
the. beam were measured at the locations shown in figure 4 with dial gages 
having a sensitivity of 0.0001 inch. In order to take advantage of this high 
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Figure 4.- Schematic diagram of test setup 
used for measuring D x . 


sensitivity, each dial gage was fitted 
with a screw with which the extension of 
the gage spindle could he adjusted until 
it just touched the surface of the beam. 
This positive screw action eliminates the 
effect of the internal friction of the 
dial gage and gives better repeatability 
of deflection measurements. Contact was 
detected by an electric circuit which 



Figure 5 .- Typical specimen and test setup used in the experimental determination of D x . 


activated an electron-ray tube. The electric detecting technique is similar 
to that described in reference 6. 


The bending stiffness Dx was calculated from the following formula obtained 
from elementary beam theory for a uniform beam subjected to a constant moment: 


D x = 


(FI) 


x s ( L 


2cw c 


( 21 ) 
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where 

F load applied at each end of beam, lb 

Z distance between load and support, in. 

x s distance from left support to any station, in. 

L distance between supports, in. 

w g deflection at station x s , in. 

c width of beam, in. 

Equation (21) was applied at three stations: x s = and The three 

values thus obtained differed from one another by less than 2 percent in any test; 
the average of the three values was taken as Dx. (This procedure is similar to 
that used in ref. 1.) 

Discussion of results .- A comparison of the experimental and theoretical 
values of for the three corrugation patterns are as follows: 

For the small vee, 

(D x ) = k-2.220 ± 200 in- lb 
' 'exp ’ 

(Dx) th = ^2,230 in- lb 

For the large vee, 

(Dx) eX p = 30,170 ± 120 in- lb 
(Dx) th = 30,280 in- lb 

For the square, 

(Dx) exp = 52,850 ± 490 in- lb 
(Dx)th = 51,730 in- lb 

As might be expected, the results of the tests for D x compare very favorably 
with the theoretical values. The theoretical value of Dx fell within the range 
of experimental values for both the small-vee and the large-vee corrugation pat- 
terns. The theoretical value of Dx for the square corrugation pattern fell 
outside the experimental range by less than 2 percent (being slightly lower than 
the lowest experimental value). 
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Experimental Verification of Formula for Dy 


Specimens . - Hie three panels used to determine the experimental values of 
Dy were 32 inches long and 6 inches wide. In this case the corrugations were 

normal to the long dimension of each panel. Each specimen was fabricated of 
Inconel X sheet which had an average thickness of about 0.0106 inch and a varia- 
tion in thickness of ±3 percent. These test specimens are shown in figure 6. 



Figure 6.- Three specimens used in vibration tests for measuring Dy. 

Test apparatus and procedure . - Vibration tests were performed on three 
panels in order to determine the experimental values of Dy. Each panel was 

hung by 8-foot-long nylon cords attached at the first-mode nodal points located 
22 percent of the panel length from each end. This arrangement approximated 
free-free boundary conditions and is shown schematically in figure 7- The panels 
were vibrated with an air shaker (ref. 7), a noncontacting inductance pickup was 
used to detect the deflections, and a recording oscillograph with a paper speed 
of 7.5 inches per second was used to obtain frequencies and maximum amplitudes. 

For each panel, several natural frequencies were obtained from which values 
of Dy could be calculated with the use of the following equation: 


D 


y 



(22) 
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where 



Figure 7.- Test setup used for 
measuring Dy. 

patterns are as follows: 


p mass per unit length of panel, 

2/2 

lb -sec / in. 

o> n natural frequency of nth mode of 
vibration of a panel, 
radians/ sec 

3 n a characteristic number from beam- 

vibration analysis (numerical 
values of 3 n were obtained 

from ref . 8) 

c width of panel, in. 

Equation (22) is the frequency equation 
for beam vibration in which the beam 
stiffness has been replaced by the 
quantity cD . 

Discussion of results . - Comparisons 
of the experimental and theoretical 
values of D for the three corrugation 


For the small vee, 

(dJ) = 8.73 to 9.36 in- lb 

' ** exp 

Nth * 8 - 53 in - ll> 

For the large vee, 

(dJ) = 11.12 to 11.69 in- lb 
\ ■*/ exp 

Nth ■ lx - 26 in - lb 

For the square, 

( Dy )exp = t0 10,85 in- lb 

Nth - ln - lb 

The experimental and theoretical values of Dy are also shown in figure 8. 


14 





As shown in this figure, the experimen- 
tal values of Dy for a particular 

panel varied with frequency. However, 
each of the values of Dy was repeat- 

able within 1^ percent . Over most of 

the frequency range, the theoretical 
values of Dy are conservative, those 

for the square corrugation being the 
most cons ervat i ve . 

Experimental Verification of 
Formula for D Yir 

Specimens . - Each specimen used for 
the evaluation of D„^ r was fabricated 

of Inconel X sheet which had an average 
thickness of about 0.0106 inch with a 
variation in thickness of ±3 percent. 

As in the tests for D x and Dy, the 

three specimens tested represented the 
small-vee, the large -vee, and the square 
corrugation patterns. The panels were 
about 18 inches square. 

As was pointed out in reference 3* when testing sandwich panels to deter- 
mine torsional stiffness, the deformation must be restricted to pure twist. 

This restriction means that the panel must deform in such a way that normals to 
the middle plane remain normal; that is, no transverse shear deformation can be 
allowed. This same requirement applies to corrugation-stiffened panels, and it 
implies the need for two restrictions on the edges of the panel: 

(1) The edges must remain straight. 

(2) Each element of the edge cross section in the xz- or yz-plane must 
rotate through the same angle (the angle of twist of the panel) . 

For restriction (1), the two edges parallel to the x-axis are kept straight 
during testing by the corrugations running parallel to the x-axis. In order to 
keep the other two edges straight, it was necessary to attach aluminum plates 
along the edges of the panel parallel to the y-axis. These edge plates were 
sufficiently thin so that their torsional stiffness was negligible in comparison 
with that of the panel. The edge plates were attached to the panels with screws 
cast into plastic extending into the ends of the corrugation cells about 3 /^ inch 
The screws were inserted through holes in the edge plates, and nuts were mounted 
on the screws to allow the plates to be tightened against the corrugation cells. 
The tightness of these nuts in some cases plays a significant role in the achieve 
ment of restriction (2), and this role is discussed in detail in a subsequent 


i2r 
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Figure 8.- Experimental and theoretical 
values of D v . 
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section entitled "Discussion of Results." 
The three panels with screws cast in the 
corrugations and aluminum edge plates 
are shown in figure 9- 

Test apparatus and procedure . - The 
panels were supported at diagonally 
opposite corners , and an aluminum rod 
was suspended from the other two cor- 
ners. Loads were applied by adding lead 
shot to a bucket hanging from the center 
of the aluminum rod. 

Nine dial gages having a sensitivity 
of 0.0001 inch were used to measure the 
twist of the panel. The electrical sys- 
tem of measuring deflections noted in the 
section entitled "Experimental Verifica- 
tion of Formula for D x " could not be 

used here because of the instability of 
the panel on its supports. When the 
electrical technique was used, small dis- 
turbances, such as dial gage contact and 

Figure 9- Specimens (including edge vibrations, caused the panel to rotate 

plates) used to measure D xy - slightly about the line connecting the 

two supports during individual dial -gage 
readings. However, when all the dial 

gages made contact with the panel at the same time, the panel was stabilized and 
did not rotate while the dial gages were being read. The forces applied by dial- 
gage contact were nevertheless negligible. Figure 10 shows a typical specimen 
and test setup used in the torsion tests. 

Loads were applied in 5-pound increments up to a maximum of 25 pounds. 
Deflections of the panel were measured at the nine points shown in figure 11. 

From the measured deflections, which were linear across the panel in the x- 

c^w 

and y-directions, the angle of twist per unit length was computed. The 

dx dy 

twisting stiffness was then obtained from the formula 




where the twisting moment M^y is equal to one-half of the torque divided by 
the width of the cross section resisting the torque, or, equivalently, 
is equal to one-half of the load acting at each corner. (See ref. 5 *) 




I 



IT 




ft 


Discussion of results . - Comparisons 
of the experimental and theoretical val- 
ues of D X y are as follows: 

For the small vee, 

(V)exp * 5 ’ 125 * 50 ln - 11J 

(V)th ' k ’ lk9 in ' lb 

For the large vee, 

(D xy ) exp = 5,837 ± 15 in-lb 


Figure 11 .- Location of dial gages on / j) \ _ c 034 in-lfc 

torsion specimens. \ Xy /th 

For the square. 


( D *y)exp ■ 5 ’ 75 ° ln ‘ lb 

K) th - in - ib 

For the small -vee corrugation pattern the theoretical result is conservative by 
7 percent, and for the large-vee pattern the difference between theory and experi- 
ment is less than 2 percent. 

A more thorough discussion is required in order to understand the results of 
the tests on the square corrugation pattern. The stiffness measurements on the 
vee-shaped corrugations were essentially independent of whether the nuts on the 
edge-plate mountings were tight or loose. For the square corrugation pattern, 
on the other hand, large variations in the apparent torsional stiffness-*- were 
measured and depended on the tightness of the nuts. For instance, in tests in 
which the edge -plate connections were loose, values of the apparent torsional 
stiffness as low as 2,900 in-lb were obtained on the square corrugation specimen. 
The experimental value reported above, 5,750 in-lb, was obtained from a specimen 
having tight edge -plate connections and having screws cast in every cell, rather 
than in every other cell as shown in figure 9- 

The square pattern falls into a category of shapes for which special pre- 
cautions must be taken to comply with restriction (2) discussed under "Specimens." 
The flat portion between the corrugation cells appears to be significant in this 


^-Since D xy is a stiffness associated with a deformation which conforms to 

restrictions (l) and (2) stated in the section entitled "Specimens," a measured 
stiffness based on overall deformations which do not conform to the restrictions 
is referred to as the "apparent" value of D^. 
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regard. Under torsional loads the defor- 
mation of the flat portion may be par- 
tially independent of the deformation of 
the cells. The panel as a whole, there- 
fore, may have a larger angle of twist 
than the individual cells. This situa- 
tion is illustrated in sketch 1 of a 
typical cross section of a panel having 
a flat portion between cells. The dashed 
outline denotes the undeformed shape and 

the solid outline denotes an exaggerated deformed shape emphasizing the possible 
difference in behavior between the flat portions and the cells. With the dial- 
gage techniques used in this investigation, it is likely that angle (l) is meas- 
ured although the cells rotate only through the smaller angle (?). 

The edge plates provided a reasonably effective means for eliminating the 
independent action just discussed. The connections between the panel and edge 
plates were such that nuts could be used to tighten the edge plates against the 
corrugation cells. When the connections were tight, rotation of the edge plates 
forced the cells having connections to rotate also, but, because of unavoidable 
slippage, not necessarily through the same angle as the edge plates. This slip- 
page, even when the connections were tight, may account for the fact that the 
experimental value of I^y for the square corrugation specimen was 9 percent 

lower than the theoretical value. 

Since the stiffness measurements on the panels with the vee-shaped corruga- 
tions were not affected by the degree of tightness of the edge plates, it 
appeared that additional experimental information on the influence of flat por- 
tions between the cells could be obtained by testing a panel consisting of 
alternating vee-shaped corrugations and flat portions. Consequently, a panel was 
tested which had a cross-section configuration as indicated in sketch 2. The 

external dimensions of this panel were 
similar to those of the other panels 
used in the torsion tests. The result 
of this test was that the apparent tor- 
sional stiffness with connections loose 
was about three-fourths of the apparent 
torsional stiffness with connections 
tight. This result substantiates the 
idea that the partially independent 
action of the flat portions and the 
cells is an important factor in the tor- 
sional behavior of corrugation-stiffened panels which have a flat portion between 
corrugation cells. 


"i * ' 
* 


' / 

\ / 

'---Rows of spot welds 


Sketch 2. 





Sketch 1. 
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CONCLUDING REMARKS 


This report presents theoretical formulas for five orthotropic elastic con- 
stants associated with pure bending and pure twist for a corrugation-stiffened 
panel. The constants are the bending stiffnesses parallel to and normal to the 
corrugation direction (d x and Dyj, the twisting stiffness ^D X y) , and the 

Poisson's ratios associated with bending ( p x and Py)- Theoretical values of 
and D are compared with experimental results for three corrugation 

y 

patterns — two vee shapes and a square shape — and generally showed good agree- 
ment. It is shown that for corrugation-stiffened panels typical of practical 
construction, the product u d.. is very small compared with unity. 

x «y 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., November 21 , 1962. 



APPENDIX A 


DERIVATION OF FORMULA FOR Dx 


The panel Lending stiffness D* is simply the Lending stiffness per unit 
width of a Learn having the cross-sectional shape shown in figure 2. When the 
panel Is considered to Le made up of a series of cells and the contributions of 
each element of the cross section are totaled, the following equation Is obtained: 


2Ep 

(a + b)t p 5 

P 

12 


2E C 

(a + i)t c 5 

P 

12 


\ J 

+ ' t c r 1^2 + + 2t c eqrq 2 sin 9 + eq 2 rq0t c + ^ k.3sin 2 0 + cos 2 0 


12 


+ ( t p + ^- + j + | sin 0 - yj kt c 


+ t c r 2 5j| + sinjej + 2t c e 2 r 2 2 sin 0 


e 2 2 r 2 0t c + + (e 2 + r 2 ^t c g 


(Al) 


where 


%-(L+a) + (tp+^)(a + i)t c + ^tp + ^ + rq- ^ ~ -j r^tc 


y = - 


n i t P rp sin 0 

+ I t p +-|. + j+^L_Jjkt c +lt p + -|.+d - r 2 + — — g lr 2 0t c + lt p +-^ + djgt c 


^ tp(b +a) + (a + i)t c + rq9t c + kt c +r 2 0t c + gt c 
Ec * ‘ *“ 


(A2) 
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APPENDIX B 


DERIVATION OF FORMULA FOR D y 


A single cell of a corrugation- stiffened panel under pure bending is shown 
in figure 12(a). As indicated, each cell must carry the total bending moment My 


applied to the panel. In this derivation, the term "plate" refers to the flat 
cover plate to which the corrugation is attached and the term "panel" refers to 
the composite structure. At the weld line, where the plate is joined to the cor- 
rugation, the moment is distributed to the corrugation and to the plate in pro- 
portion to their respective stiffnesses. The plate and the corrugation are now 
treated as separate structures each having its own force system and deformation 
as shown in figure 12(b). The curvature in the x-direction resulting from My 



( a) Undef ormed 
composite cell 
with My loading. 



(b) Loads and 
deformations of 
cell components. 


is considered to be zero. The forces 
and deformations are related in the 
following way: 

(1) The slopes of the plate and 
corrugation are equal at the weld line. 

(2) - The extensions of the plate 
and corrugation are equal at the weld 
line . 


Figure 12.- Loading on corrugation ceil. ( 3 ) ^ moment carried by the cor- 

rugation plus the moment carried by the 
plate equals the moment carried by the panel; that is. 


M c + Mp = My (Bl) 

(4) The slope at the weld line for each element is the same as the slope of 
the panel which is determined by Dy and 

(5) Equal and opposite (in-plane) loads P act at the weld line. Symmetry 
and overall equilibrium of the plate or corrugation normal to the panel dictate 
that no normal forces act at the weld line. 



The corrugation is analyzed first. The bending strain energy U caused by 
the loads P and moments M c is computed. The slopes and extensions caused by 
bending are then obtained by taking the appropriate partial derivatives of the 
strain energy. The strain energy is calculated from 



/ 


(Moment) ds 


(B2) 


where D c is the plate bending stiffness of the corrugation material, and the 
integration is taken over the corrugation cross section. In equation (B2), 
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rather than (El) is used because the bending is considered to be cylin- 


drical, that is. 


a 2 w 

dx 2 


= 0. If the corrugation is divided into elements AB, BC, 


CD, DE, and EF, as shown in figure 13, the contribution of each corrugation 

Element to the strain energy is given 


Figure 13.- One-half of a single-weld- 
line-type corrugation cell 
(appendix B) . 



by the following equations for the case 
in which the panel has only a single 
weld (a = 0) : 


For element AB, 


U = So 


2 

M dy 
c 


(B3) 


For element BC, 
, 8 


U = 2D~ J 0 fc-^ 1 " 008 ®^J 2r i 

(B4) 


For element CD, 


k 2 

0 - afc J 0 [ H = - P(J + s sin 9 >1 ds 


(B5) 


For element DE, 


U = 





- 9) 



For element EF, 



(B6) 


(bt) 
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sin 20 


Expanding, integrating, and totaling these contributions gives 
U = ^j^M c 2 i + Tl M c 2 0 - 2M c Pr 1 2 (0 - sin 0) + P 2 ^ 5 ^- - 2 sin 0 

+ Mc 2 k - 2PMc^k + sin oj + P 2 ^^ + jk 2 sin 0 + ^ sin 2 ©j + r 2 M c 2 0 

- 2r 2 M c P^n0 - r 2 0 cos 0 + r 2 sin 0) + r 2 P 2 |n 2 0 - 2nr 2 0 cos 0 + 2nr 2 sin 0 
+ (39 - 3 sin 0 cos 0-20 sin 2 ©) + M^g - 2M c Pgd + P 2 gd 2 | (B8) 


According to Castigliano ' s theorem, the slope of the corrugation at the weld 


line is given by - 


dU 

SMc 


where the sign is determined from the coordinate system 


shown in figure 13* The slope of the plate at the weld line is given by - 


Mpb 

Dp 

and, because of the weld, is equal to the slope of the corrugation at the weld 
line. Thus, 


Dp c3Mc 


= Mci + r^Mc© - Pr x 2 (0 - sin 0) + l^k - P ^jk + sin ^ 


+ r 2 M c 0 - r 2 P^n0 - r 2 0 cos 0 + r 2 sin 0) + ^g - Pgd 


(B9) 


Another equation can be obtained by equating the extensions of the corruga- 
tion and plate at the weld line. In this analysis, only the corrugation exten- 
sion caused by bending is considered and the plate extension is neglected. This 
assumption seems justified because extensions caused by direct strains in the 
plate are of the same order of magnitude as the extensions caused by direct 
strains in the elements of the corrugation. The direct strains in the corruga- 
tion elements are negligible compared with the bending deformations. However, 
when the depth d of the corrugation is small, the extensions caused by direct 
strains have to be considered. 


2k 



T7- 


The extension of the corrugation at the weld line is given by — . 

dP 

Equating the extensions of the corrugation and plate at the weld line gives 


dU _ _1 

dP D c 


-r.^0 - sin 0) + Pr^ 5 ^- - 2 sin 0 + sl ? 29 j - M c ^jk + sin 0^ 
+ P ^) 2 k + jk^sin 0 + ^- sin 2 0^ - ^n0 - r 2 0 cos 0 + r 2 sin 0^ 


+ r 2 p|n 2 0 - 2nr 2 0 cos 0 + 2nr 2 sin 0 + ^30 - 3 sin 0 cos 0-20 sin 2 0^ 


- M_gd + Pgd 


= 0 


(BIO) 


Grouping the coefficients of P and M c yields 

P = HMc (Bll) 

where 


H = — 


r l 2 (0 - sin 6 ) + jk + sin 0 + r 2 (n0 - r 2 0 cos 0 + r 2 sin 9^ + gd 


Tt^(— - 2 sin 0 + — 2 ^ \ + j^k + jk^sin 0 + —■ sin^0 + r 2 n^0 


L 1 2 


2 nr -3 9 cos 9 + 2m>> sin 0 + 


r 2' 


■(30 - 3 si 


sin 0 cos 0-20 sin^0 


I'j + gd^ 


Substituting HMc for B in equation (B9) gives 


Me Dp 

"P-V5?* 


(B12) 


(B13) 


where 


X = i + 9^r^ + r2 j + k + g - H [jl 2 ( 0 - sin 9) + jk 
+ sin 0 + r 2 ^n0 - r 2 0 cos 0 + Tq sin 0 ) + gdj 


(Bl4) 
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If the expression given in equation (B13) is substituted for Mp in the 
equation 


the result is 


M c + Mp = My 



(B15) 


Because the curvature in the x-direction is considered to be zero, the moment- 
curvature relationship for the panel is 


- 


Py a 2 w 
1 - M-x^y dy^ 


(Bl6) 


By considering to be constant and by substituting the expression for 

«y given by equation (B15) into equation ( Bl6 ) , the slope of the panel at the 
weld line is given by 



(B17) 


Equating the slopes of the corrugation and panel at the weld line gives 


Me 


X 




(BIB) 


where Dy' is the bending stiffness per unit 
sion a equals zero, that is for the case in 

line instead of a double weld line. Equation 
to give 


width of the panel If the dimen- 
which the panel has a single weld 

/ B ' 

(Bio) can be solved for - — & 

1 - HxHy 


Dv' bB c 

— — = + Bn 

1 - JJ^ly X TP 


(B19) 


For the case in which the panel has a double weld line, the dimension a 
does not equal zero, and the segment of length 2a has to be considered. This 
segment is treated as a composite beam. The location of the neutral axis for 
such a beam Is as follows (see ref. 9 ): 
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(B20) 



where z is the distance from the outer surface of the plate to the neutral 
axis as shown in sketch 3 . 


Neutral 



Sketch 3* 


The bending strain €y at a distance z 
from the neutral axis is given by 



(B21) 


where r is the radius of curvature in 
«/ 

the y-direction. 


Again, when bending is considered to be cylindrical, the bending strain e x 
equals zero. Therefore, with the use of Hooke's law and equation (B2l), the 
bending stress in the y-direction is given by 


CT,r = 


Ez 

(l - p 2 )r y 


The normal stresses distributed over the ends of the segment must produce 
couples equal to the applied moment, or 


(B22) 


% ■ 

J 5-tp-tc 


cr c z dz + 



CTpZ dz 


(B23) 


where xj c is the stress in the y-direction acting on the corrugation and Op is 
the stress in the y-direction acting on the plate. 


Substituting equation (B22) into equation (B23) and using the appropriate 
values of E and p. gives 


*v - 



L 


z-tp-t c 


z 2 dz 




(B24) 


By comparing equation (B24) with the well-known moment- curvature relationship 
for cylindrical, bending of isotropic plates, that is 


“sr- 



(B25) 
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the tending stiffness of the segment between weld lines can be defined as 


r *» Z 2 dz + r ^ 

1 " ^xM-y 1 - (i c 2 « z-tp- t c 1 - (ip 2 z-tp 


(B26) 


where Dy" is the bending stiffness per unit width of a beam cut from the seg- 


ment of length 2a. After performing the integration. 


T) n 

D y 


1 - dxdy 


is given as 


v 


E r 


1 - (1-jjPy 3^1 _ H C 2 

Ep 




y ^c^ 3tpt c + 3tp t c + 3z 

) 


■“ 6ztp"tQ “ 3^^ 




t p 3 - 3zt p 2 + 35 2 t p 


(B27) 


For the case in which the cover-plate material is the same as the corruga- 
tion material. 


D y" _ Sfrc + t p’)^ 
1 - dxPy 12(l - (i 2 ) 


(B28) 


The change in slope between the end and the center of the segment of length 2a 
is given by 



(B29) 


The change in slope between the weld line and the center line of the corrugation 
cell is 


Myb 



(B30) 


The total change in slope between the center of the segment of length 2a and the 
center line of the corrugation cell is the sum of expressions (B29) and (B30): 



(B31) 
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Dividing expression (BJl) by one-half the pitch gives the change in slope 
per unit length. Since the change in slope per unit length is also equal to 


My(l ~ ^y) 

Dy 


there results 



The bending stiffness D y is then given by 


(B52) 


Dy = _/ ■ + V (B55) 

V D y' 

where Dy" is given by equation (B27) or (B28) and D y ’ is given by equa- 
tion (B19) . It is shown in appendix D that for a corrugation- stiffened panel, 
the product M-xP-y I s usually small compared with unity and can be neglected in 

the computation of Dy' and Dy" . 
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APPENDIX C 


DERIVATION OF FORMULA FOR D xy 


It was mentioned in the body of this report that a corrugation- stiffened 
panel might be considered as a series of cells. Each of these cells has the 
torsional rigidity GJ of a thin- walled tube. The torsion- twist relationship 
for a thin- walled tube is given by 


T = GJ9 Z 


(Cl) 


where 

T total torque on a cell 

GJ total torsional rigidity of a cell 

angle of twist per unit length of a cell 

The torsion- twist relationship can also be expressed as 

T = GJ0 i (C2) 

Where the bars indicate that the quantities T and GJ are divided by the width 
of the cross section considered. 

The torsion- twist relationship for an orthotropic panel is 

M X y = Dxy (C3) 

The torque per unit width T is considered to act on only two edges of a rectan- 
gular orthotropic panel whereas, to be consistent with the differential equation 
for lateral deflections, the twisting moment M xy is considered to act on all 
four edges. Therefore, for the same deformation to take place, the torque per 
unit width T must equal twice the twisting moment M xy . A comparison of equa- 
tions (C2) and (C3) then yields 

Dxy = f (04) 


where 



Torsion constant per unit width of cross section (from Bredt 
theory, see ref. 10) 


(C5) 


30 



in which 


A area enclosed hy corrugation cell 
(see sketch 4) 

ds elemental length along center line 
of sheet making up corrugation 
cell 

t thickness of sheet associated 
with ds 

p pitch of corrugation pattern 



Sketch 4 . 


The integral is taken along the center line of the corrugation cell. The 
shear modulus G is retained inside the integral to account for the possibility 
of different materials for the plate and corrugation. The area A is given by 


A = b^t c + tp^ + d jb + g - i + ^r 2 - r^ tan |j + 2 (Vj_ 2 _ r 2 2 Jtan | + ^r 2 2 - r^ 0 

(C 6) 

The integral along the center line of the corrugation cell is given by 

§ m * sp; + + r i e + r 2 e + k + 8 ) (C7) 


Substituting equations (C5) to (C7) into equation (c4) gives 


D - 3L 
^xy - ~ 


^•|b(tc+tp) +djb + g-i + (r 2 -r 1 )tan |j + 2 ^2 - r 2 2 j tan § + (r 2 2 - r x 2 ) ej 


+ t^( i + r l 0 + r 2 e + k+s ) 


(C8) 
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APPENDIX D 


DERIVATION OF FORMULAS FOR p x AND n y 


For small deformations, p x may be defined as the ratio of the panel curva- 
ture in the y-direction to the panel curvature in the x-direction when the moment 
M x is the only loading on the panel. The definition of p y is obtained by 
interchanging x and y in the definition of p x . Just as in the case of iso- 
tropic plates, these definitions of n x and p y are applicable only when the 
deformations are small. 


An expression for p x is obtained by use of a method which has some similar- 
ities to that used in the derivation of Dy. Consider a rectangular corrugation- 
stiffened panel and assume the plane x = 0 passes through the center of the 
panel. (See fig. 1.) A moment M x is assumed to be applied to the panel and 
yields a curvature K in the x-direction. The curvature in the yz-plane, which 
is independent of x, is calculated and this curvature divided by K is the 
Poisson's ratio p x . In order to calculate the curvature in the yz-plane due to 


M x , the corrugation and plate which comprise the panel are first considered to be 
separated at the weld lines and to act independently. The axial- strain distribu- 
tion caused by M x is applied to the corrugation and plate, and their deforma- 
tions are calculated at x = 0. It is found that redundant forces and moments 
are necessary to match up the extensions in the y-direction and the slopes in the 
yz-plane at the weld lines. These forces and moments are also determined by the 
use of Castigliano ' s theorem as was done in the derivation of D y . The other 


bending Poisson's ratio p. y can be determined from the relation 
can be derived from the theorem of reciprocity. (See ref. 3-) 


Py = Px 
D y D x 


which 


Deformation of Corrugation 

A moment M x acting on the corrugation, figure l4(a), is assumed to give a 
linear strain distribution e x as shown in figure 14(b). If the corrugation 
and plate are separated, moments M and loads P act at the weld line as shown 
in figure l4(c). The corrugation, which is considered to have a single weld line 
at present, is analyzed first. Figure 15 gives the dimensions of the corrugation 
cell and divides the cell into elements AB, BC, CD, DE, and EF. For all elements, 
the axial strain e x is considered to be uniform through the thickness t c or 
tp and equal to the axial strain of the middle surface of the element. The 
deformation of the corrugation caused by the axial strain e x is now considered. 

Element AB .- The axial strain e x in element AB is given by 

e x “ " K ( e l + r l) ( D1 ) 
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** - Kz 


I 


z 


'-■Neutral axis 



Figure 15 . - One-half of a single-weld- 
line-type corrugation cell 
(appendix D) . 


where K is the primary curvature 
resulting from Considering cTy 

to he zero and using the stress-strain 
relations for plane stress gives 

€ y = K ^c( e l + r l) (D2) 



(a) Undeformed composite 
cell with M x loading 

(b) Strain distribution 
resulting from M x . 

(c) Loading at weld line 
and deformations of 
cell components. 


The extension of element AB in the 
y-direction is then given by 

= Kip^e! + r^ (Dj) 

Element BC .- Element BC is analyzed 
as a sector of a cylinder deformed by a 
uniform axial load and a bending moment. 
Coordinates and displacements are shown 
in sketch 5. 



Figure 14.- Loading on corrugation cell. 


Sketch 5. 



The strain- displacement relationships for element BC are 


e 


7 


e =Su 
x ax 

(D4) 

_ 1 Sv w r 
r l c >0 " r l 

(D5) 

_ l 5u , dv 
r l § 0 " dx 

(D 6 ) 


The strain e x is given by 

e x = ~ K ( e i + r i cos 0) (D7) 

Thus, from equation (d4) 

u = -Kx^ej + ^ cos 0^ + f( 0 ) (d8) 

At x = 0, u = 0; therefore, f(0) = 0. 

With the assumption of Ja = 0, equations (d 6 ) and (D 8 ) give 


v = v x ( 0 ) - ^ sin 0 


(D9) 


where v^_(0) is an arbitrary function of 0 only. Since all strains are 

assumed to be functions of 0 only, examination of equations (D5) and (D9) leads 
to the conclusion that v r must be of the form 


w r = w l( 0 ) 



(DIO) 


where w^(0) is an arbitrary function of 0 only. Since a^ = M^ = 0 at 0 = 0 
and 0 = 0 , it can be shown from the equations of equilibrium that CT 0 = = 0 

everywhere in this element. This information provides two relations from which 
expressions for v^(0) and w-^(0) can be obtained. By using = 0 and the 

stress- strain relations for plane stress, there is obtained 

e 0 + ^c e x = 0 (Dll) 
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and substituting for e0 and e x from equations (D5) and (D7)> respectively, 
yields 


i ^ - ?r - ^ K ( e i + n oos *) - 0 


(Dia) 


Substituting the formulas for v and w r from equations (D9) and (DIO) , respec- 
tively, into equation (D12) gives the first relation: 


7? ^ - ^ K ( e l + r l CQS ^ = 0 


(D13) 


The second relation is obtained from M0 -= o and from the moment- curvature 
relation 


where 


M0 - -Dc^0 + - 0 


y . _ 1 fdv . Sfwi 

$ a0 2 


X Y = 


a 2 


w r 


dx 2 


(di4) 


(D15) 


(di6) 


Substituting equations (D15) and (Dl6) into equation (Dl4) gives 




1 fdv . d v r | _ 

r^W 50 2 / ° dx 2 


(D17) 


Again using equations (D9) and (DIO) for v and w r , respectively, yields the 
second relation: 


“rfljr " *2“ C ° S ^ + S?" + COS = ^c* cos 0 


B0 2 


(DIB) 


dv n 


Substituting the formula for — from equation (D13) into equation (Dl8) gives 

c^w-. 


wi + (JL c riK^e]_ + r^ cos 0) + 


d0 c 


= ^ C K cos 0 


(D19) 


35 



or 


d^wq 

S0 2 


W 1 


-^c r l Ke l 


(D20) 


When equation (D20) is solved for Wq and the result is substituted into equa- 
tion (DIO), the following equation results: 


w T 


= cq sin 0 + C 2 cos 0 - Hc r l Ke l - cos 0 


(D21) 


Similarly, using the solution to equation (D20) in conjunction with equa- 
tions (D13) and (D9) yields 


v = -cq cos 0 + C 2 sin 0 + (i c Krq^sin 0 - 


Kx3 


sin 


0 + c- 


(D22) 


where cq, cp* and cj are arbitrary constants representing rigid-body dis- 
placements in the x = 0 plane. The quantities of interest in this analysis are 
the relative rotation of one end of element BC with respect to the other end in 
the x = 0 plane and the extension of the element in the y-direction in the 
x = 0 plane. These quantities can be obtained by ignoring rigid-body displace- 
ments. By setting Cq = C 2 = c-z, = x = 0 in equations (D2l) and (D22), it can be 

seen that there is no relative rotation of the ends of the element; that is, 

= 0. The extension of the element in the y-direction is 

x=0 


dw r 


52 - -(sin 0)w I 


0=0 

x=0 


+ (cos 9)v 


0=9 

x=0 


(D23) 


or 

p 

6q = q c rqKeq sin 9 + p c Krq sin 9 cos 9 

Element CD .- The axial strain in element CD is given by 

e x = Kz 

By considering the in-plane stress in the direction to be zero 
using the stress- strain relations for plane stress, there results 


{D2k) 

(D25) 
o) and by 


eg = -Hc Kz 


(D26) 
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The total elongation In the |-direction from point C to point D is given by 


5 t = / c ” n d| 

The extension in the y-direction is given by 

= 5| eos 6 


(D27) 


(D28) 


or 



H C K 


cot 9 


/ \2 

/ \2 

( e l + r l cos 9 j 

- ( e 2 + r 2 cos 9 j 


(D29) 


Element DE .- The equations for v and w r for the lower arc (element DE) 
are similar to those of the upper arc (element BC). The only differences are the 
substitutions of -e 2 for e^ and r 2 for r^_. Angle 0 is measured from the 

same reference line as for the case of the upper arc. Displacement w r is posi- 
tive toward the center of curvature of element DE. When rigid-body displacements 
are ignored. 


w r = Pc r 2 Ke 2 


Kx £ 


cos 0 


v = p c KT2^sin 0 - sin 0 


(D50) 

(D31) 


As before, the relative rotations of the ends of the element are zero. The 
extension in the y-direction is given by 


5]^ = -(sin 6)w r 


0=9+11 

x=0 


+ 


(cos e)v 


0=9 +TT 

x=0 


(D32) 


or 


O 

6]^. = -p c r 2 Ke 2 sin 9 - H c Kr 2 sin 9 cos 9 


(D33) 


Element EF.- The axial strain for element EF is given by 


6 x = K ( e 2 + r 2) 

If 0 y is considered to be zero, 

e y = -^c K ( e 2 + r 2j 


(D3^) 


(D35) 
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The extension 6^ of element EF is then given by 

S 5 = -^ c Kg(e 2 + r 2 ) 


(D5 6) 


Analysis of Redundant Loads and Moments 


In order for the corrugation and plate to have equal slopes and extensions 
at the weld line, equal and opposite redundant moments M and loads P act at 
the weld line. The deformations caused by these loadings are considered to be 
independent of the deformations caused by the linear strain distribution e x . 

The extension of the corrugation in the y- direction caused by moments M and 


loads P 


is given by 


where 

Bp 


u 


is the strain energy of the corrugation due 


to bending. 


An expression for 


oU 

BP 


is given by equation (BIO). 


The total extension of the corrugation & c is given by 

5 




(D37) 


n=l 


The extension of the plate Sp is given by 

5p = (XpKb^y - (D58) 

Equating the extensions of the corrugation and plate at the weld line gives 

6 

&n + ^ = 0 

^ Bp 


(D39) 


n=l 


where 


5 6 = - 5 P = -^p Wy - 


\ 

i n equation (BIO), Mq is replaced by M then can be expressed 

Bp 

linear combination of M and P as follows: 


(d4o) 


as a 


I “ 5^™ + 4P > 


(d4i) 
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where -A is the numerator of H and $ is the denominator of H. (See 

ecL- (B12).) 


Let 


o 

X 8n = 


KE 


(D^2) 


n=l 


where 


E = iia^ej^ + r-[J - gp c (e2 + r£) + P c r l e l si h 9 + M- C r l 2sin 9 cos 9 
cot 0 [7 \2 / \?1 

+ — — (e p + r p cos 0 J " ( e 2 + r 2 cos 9 ) “ • 1 c r 2 e 2 SILn 9 

- (i c r 2 2 sin 0 cos 0 - p p b/V - ^ 


(Dl+5) 


Upon substitution of equations (D4l) and (d 42) into equation (D39), the following 
equation is obtained: 


KE = -=i-{AM + $>P) 
D c 


(d44) 


or 


D C KE + AM 

<& 


(1*5) 


Since ^ = -H, equation (d 45) becomes 


P = + HM 

$ 


(d46) 


Inasmuch as the deformation caused by the linear strain distribution due to 
M x makes no contribution to the slope, the slope of the corrugation at the weld 

line (relative to the slope at the center line of the cell) is given by 

oM 

where a counterclockwise rotation is considered positive. An expression for 

^ is given by equation (B9) by replacing by M. The slope of the plate 

oM 

at the weld line is given by Equating the slopes of the corrugation and 

bp 

plate at the weld line gives 


c)U = Mb 

SM “Dp 


(1*7) 
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From equation (B9), 
as follows : 



can be expressed as a linear combination of 


M 


and 


P 


^ = -^(QM + AP) 
SM Dc 


(d48) 


where 


fl = i + 9(r 1 + r 2 ) + k + g 


(D49) 


Substituting the formulas for P and ^ from equation (D46) and equa- 

oM 

tion (D47), respectively, into equation (D 48) yields 


Solving equation (D50) 



/ DpKZ \ 

HM + 

A (" « + ^Jj 

M 

XDp 

gives 

M 

A z/$ 

KBp 

+ AH) + 
u c 


Mb 

Bp 


(D50) 


(D51) 


Equation (D5l) can be simplified by noting that ^ = -H and from equation (Bl4) 

Q, + AH = X. Since is the ratio of anticlastic to primary curvature, equa- 

KDp 

tion (D51) determines p x if the panel is of the single- weld- line type. For the 

purpose of distinguishing this type from a double- weld- line type, the Poisson's 
ratio given by equation (D5l) is designated p x lJ "that is. 


^xl = 


-HZ 

Bp 

— X + b 

D c 


(D52) 


where H, Z, and X are given by equations (B12), (d 43), and (Bl4), 
respectively. 

If the panel is constructed with a double weld line, the segment between 
weld lines must be considered. (See fig. 2.) The extension of one-half the 
cover plate between weld lines is given by 


/ t p + t c l 

PpKa e i + r-L + , 


(B53) 
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The extension of one-half the corrugation between weld lines is given by 


licKa^ + rp) (B5 1 *-) 

The curvature of this section is simply the difference between the extensions of 
the plate and the corrugation divided by the distance between the center lines 
of the plate and corrugation and by the distance a. The ratio of this curva- 
ture to K is designated p x 2 and is 


Hx2 


- ^p 


+ 2 (^P - M(fl + rp) 
"tp + t c 


(D55) 


The ratio of anticlastic curvature to primary curvature for the entire panel is 
then 


^x = 


b ^xl + a ^x2 
b + a 


(D 56 ) 


where p^ and are given by equations (D 52 ) and (D55)> and a and b are 

dimensions shown in figure 2. By using the reciprocity theorem for elastic 
structures, p y is obtained as follows: 

Py = P x ( D 57) 

For a corrugation- stiffened panel, the ratio of Dy to D x is usually very 
small so that py is very small. Although p x may not be small compared with 
unity, the product p x py is usually small compared with unity. For example, the 
product PxM-y is less than 5 X 10”^ for the test panels considered in this 
report (see table II ) . 
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